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Abstract. We analytically study the ground-state quantum phase transition between the Haldane phase 
and the large-D (LD) phase of the 5 = 1/2 ferromagnetic-antiferromagnetic alternating Heisenberg chain 
with on-site anisotropy. We transform this model into a generalized version of the alternating antiferromag- 
netic Heisenberg model with anisotropy. In the transformed model, the competition between the transverse 
and longitudinal bond alternations yields the Haldane-LD transition. Using the bosonization method, we 
show that the critical exponents vary continuously on the Haldane-LD boundary. Our scaling relations 
between critical exponents very well explains the numerical results by Hida. 
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§1. Introduction 

Since Haldanc's prcdiction[l,2], the integer spin chains have been attracting much attention. Hida 
[3,4] tried to elucidate the properties of the spin-1 Heisenberg chain from the standpoint of the spin-1/2 
ferromagnetic-antiferromagnetic alternating chain. He considered the Hamiltonian [4] 

N N N 

nD = 2jJ2 S2i • S2,-+i + 2 J' J2 S2i-i • S2,- + Dj2{Sij-i + Sljf , (1.1) 

i=i i=i j=i 

where is the spin-1/2 operator. We assume that J > (antiferromagnetic) and J' < (ferromagnetic) 
and we are interested in only the ground-state. Hereafter we call the model (1.1) "D-model" following Hida 
[5]. 

In the case of J' — oo, where the spins 82^-1 and 82^ form a triplet, the model (1.1) reduces to a 
spin-1 antiferromagnetic Heisenberg chain with on-site anisotropy: 

J N N 

nt' = 2 E Si • Si+i + D Y^is^r , (1.2) 

3 = 1 3 = 1 

where Sj = 82^-1 + 82^ is the spin-1 operator. When J' = £) = 0, on the other hand, Hamiltonian (1.1) 
reduces to the two-spin problem. Its ground-state is the complete dimer state in which the spins 82j and 
82^-1-1 form a singlet pair (l/-\/2) (t2ii2i-i-i — i2iT2i-i-i)- When D = 0, therefore, Hamiltonian (1.1) smoothly 
connects the dimer state of the spin-1/2 chain and the Haldane state of the spin-1 chain. 

Hida [3] numerically diagonalized the finite system of Hamiltonian (1.1) in the D = case to find that 
there is no evidence for the phase transition of the ground-sate between J' = and J' = —00. Therefore 
he concluded that the Haldane phase of the spin-1 chain can be interpreted as the special case of the dimer 
phase of the spin-1/2 chain. This picture was supported by the work of Kohmoto and Tasaki [6], and Takada 
[7] who used the nonlocal unitary transformation. Hida [4] also performed the numerical diagonalization for 
the D ^ case to draw the phase diagram of the £)-model on the J' — D plane. He found the Haldane phase, 
the large-D (LD) phase, the Neel I phase and the Neel II phase. The Neel I state is the usual Neel state 
of the spin-1/2 chain, whereas the Neel II state is Hke the | • • • T2j-iT2ji2j-i-ii2j-i-2 • • •) state which becomes 
the usual Neel state of the spin-1 chain in the limit of J' — > —00. Hida [8] also numerically estimated the 
critical exponents the Haldane-LD boundary. 

Yamanaka, Hatsugai and Kohmoto [9] investigated the model 

N 

Hx = 2J ^^(5'2j<S'2j_|_i + 5'2j5'2j_|_i + AS'|jS'|j_|_i) 
3=1 

N 

+ 2J'^ 82^-1 • S2j , 

where the parameter A represents the interaction anisotropy of the antiferromagnetic bonds. We call this 
model "A-model". This model is equivalent to the spin-1 XX Z model 

I ^ 

^A=' = 2 E(^^%i + ^Pj+i + ^^J-^+i) ' (1-4) 

in the limit of J' — > —00. We note that the ferromagnetic bonds should be isotropic so that 82^ and 
S2J+1 form a triplet pair when J' —00. They mapped this model onto the highly anisotropic version of 
the two-dimensional Ashkin- Teller model and performed the "high temperature expansion" to obtain the 
ground-state phase diagram. Their phase diagram consists of several phases, including the Haldane phase, 
the Neel phase and the XY phase. Hida [5] also drew the phase diagram of the A-model by mapping the 
A-model onto the D-model. His conclusion was consistent with that of Yamanaka et al.[9]. 

Okamoto, Nishino and Saika [10] have shown that the Haldane-LD boundary of Hida's phase diagram 
(obtained by the numerical calculations) for the £>-model can be semi-quantitatively reproduced by an 
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(1.3) 



analytical method. In this paper, we proceed their work to investigate the; critic;al properties of the Haldane- 
LD transition. We show that the Haldane-LD transition is of the Gaussian universality class and the critical 
exponents vary continuously on the boundary. In §2, we explain the nature of the Haldane-LD transition. In 
§3 the bosonization approach to the Haldane-LD transition given and the critical properties are discussed. 
The last section §4 is devoted to discussion. 

§2. Nature of the Haldane-LD transition 

In the following we investigate Hamiltonian (1.1) assuming J' < 0. Performing the spin rotation around 

the 2:-axis, 

% = -5J, S^ = -S^, S^ = S^^ fori = 4Z,4Z + l , 

SJ = SJ , S] = S] , 5| = 5| , otherwise , ^ ' ' 

we can transform Hamiltonian (1.1) into an antiferromagnetic form [4] 

2N 

n = 2Jo Y^iiS^SU + ^J^jVi + A5|5|+r) 



where 



5± 

and 



. . (2.3) 

J'l , i_|j'|_(25_2|J'|) 



1 + |J'| 1 + |J'| 



J' EE J' /J , D = D/J . (2.4) 



Hamiltonian (2.3) can be interpreted as a generalized version of the alternating antiferromagnetic Heisenberg 
chain with anisotropy (see Fig.l). The point {\J'\,D) = (1,2) (called S point) is the solvable point where 
Hamiltonian (2.2) represents the imiform and isotropic antiferromagnetic chain. The S point is the tricritical 
point of the Haldane phase, LD phase and the Neel phase in Hida's phase diagram [4]. 

In the bosonization theory of the bond alternation problem, it has long been considered that the 6^ 
term plays an irrelevant role. However, we have recently pointed out that the ground-state of Hamiltonian 
(2.2) in the case of 5_l = and 5^ > is the dimer state [11,12]. This shows the importance of the term, 
because the ground-state would be the spin-fluid state if the Sz term played an irrelevant role. When A = 1 
and S± = Sz = I, the groimd-state of Hamiltonian (2.3) is the complete dimer state (an array of local singlet 
dimers) , and both the 6z term and the d± term give the same contribution to the excitation gap apart from 
a factor 2. Then, in A = 1 and S± = Sz > case where the ground-state is the dimer state, we expect that 
the excitation gap should behaves as 

cgap - {25^ + 5zY , (2.5) 

due to the SU(2) symmetry. Here v is the exponent of the correlation length and we do not enter into the 
logarithmic correction problem [13]. Thus we can consider that the term and the 5z term are mutually 
cooperative when 5^_5z > 0, and competing when 5±8z < 0. Of course, when A ^ 1, equation (2.5) itself 
may be no longer valid. 

Let us consider the A < 1 case, where we can exclude the possibility of the Neel state. When > 
and 5z > 0, an efFeetive singicit dimer is formed by spins and S2j+i. Their coupling is antiferromagnetic 
in the original spin (S) representation before the transformation (2.1). Since this singlet dimer is still a 
singlet dimer in the S representation, this state is the Haldane state. On the other hand, when S± < 
and Sz < 0, spins 82^-1 and S2j (their coupling is ferromagnetic in the S representation) form an effective 
singlet pair. By the transformation (2.1), this singlet dimer is transformed into the triplet dimer (like 
(l/-\/2)(t2j-ii2i + i2j-il2j) state) in the S representation. This is the = state of the spin with 5=1. 
Therefore this state is nothing but the LD state. We note that, in the S representation, both the Haldane 
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state and the LD state are the singlet dimer state, but the dimer configuration is different by one spin site 
between these two states. 

What happens in the region where < and Sz > 7 If the effect of 6z {6±) is predominant, 
the ground-state may be the Haldane (LD) state. Therefore the ground-state phase transition between the 
Haldane state and the LD state can be observed in this region. The Haldane-LD phase boundary may be 
determined by the Hne on which the effects of S± and 6z cancel out with each other. Prom equation (2.5), 
the simplest estimation for the Haldane-LD boundary may be 

2S± + 6z = 0, (2.6) 

which results in 

D = 2\J'\. (2.7) 

Although this estimation is too rough because equation (2.5) may be valid only when A = 1, the nature of 
the Haldane-LD transition is well explained by the above mentioned picture. More elaborate estimation will 
be given in the next section. 
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§3. Critical properties of the Haldane-LD transition 

Through a careful bosonization procedure, we can map Hamiltonian (2.2) onto a generalized version of 
the sine-Gordon Hamiltonian 

Wb = 2 Jo j dx{A{Vef + CP'^ + Bi cos 261 - B± cose- B^{Vef cos 61} , (3.1) 

where the commutation relation 

[e{x),P{x')]=i5{x-x') , (3.2) 
holds. The coefficients B±_ , B^ and Si are obtained directly from the bosonization procedure as 

S. = ^, 5. = ^, Bi=A, (3.3) 
a n 2a 

where a is the spin spacing. These expressions are considered to be valid when J_l-c1, ^^-Cl, A-Cl. 
Similar expressions were already obtained by Nakano and Fukuyama [14], but there was an error in the 

sign in their expressions. In their expressions, the B± term and the Bz term are mutually competing when 
5±5z > 0. From the discussion of §2, these terms should be mutually cooperative when S±5z > 0, as is 
realized in (3.1). 

Since the term Bi cos 26 is irrelevant for the A < 1 case with which we are concerned, we may neglect 
this term by setting 

Bi = . (3.4) 
For the coefficients A and C, the bosonization procedure leads to 

Of course, these expressions may be valid for A <C 1. Therefore we cannot use such expressions, because we 
require expressions valid near A = 1 to discuss the Haldane-LD transition near the S point (A = 1, S± = 
Sz = 0). When S± = 6z = 0, Hamiltonian (2.2) represents the uniform XX Z chain for which the exact 
results are available. In this case, the spin wave velocity v [15], and the power decay exponent r] [16], defined 
by i-lYiSSS^) - r-", are 



V Travel — A2 
2Jo ~ 2 cos- 1 A 



(3.6) 



2 

l + (7r/2)sin-iA ' ^^'^^ 
respectively. If we use the bosonized Hamiltonian (3.1) with Bi = B± ^ Bz — 0, we obtain 

Prom equations (3.6)-(3.8) we can immediately write down the expressions for A and C. This procedure for 
the adjustment of the coefficients was first proposed by Cross and Fisher [17] and applied by Nakano and 
Fukuyama [14,18]. If we expand A and C with respect to e = 1 — A, we obtain 
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up to the lowest order of c. Expressions (3.9) and (3.10) wore also obtained by Inagaki and Fukuyama [19]. 
Thus our bosonized expression for the spin Hamiltonian (2.2) is equation (3.1) with equations (3.3) and 
(3.10). We note that this is valid for |<5_l| < 1, < 1 and e < 1. 

To discuss the critical properties of the Haldane-LD transition, we apply the self-consistent harmonic 
approximation (SCHA), which is essentially the variational method, to the bosonized Hamiltonian (3.1) 
when < and 5z > 0. In the Haldane state, as is discussed in §2, the effect of Sz is predominant and the 
average value of with respect to the ground-state, (6), is 

(6') = , (Haldane state), (3.11) 

so that the Bz term gains the energy. In the LD state, on the other hand, the effect of Sz is predominant 
and 

{0)=w , (LD state). (3.12) 

Then we set the SCHA Hamiltonian 

Hs = 2Jo j dx{A(V0)2 + + B(j)^} , (3.13) 



with 

6 , (Haldane state) , , 

e-TT, (LD state) ' ^ ' 



where B is the variational parameter. The parameter B should be determined so that (Wb)s is minimized, 
i.e., 



dB 



= , (3.15) 



where (• • •)s denotes the average with respect to the ground-state of Tis- 

The excitation spectrum of Ws is 

ws(9) = v^^+^, , (3.16) 

where 

ql EE B/A . (3.17) 
and vqc is the excitation gap. Since Tis is harmonic, the relations 

(exp(iM</.))s =exp (^-y(02)s^ , (3.18) 

((V0)2cos(/.)s = ((V0)2)s(cos<^)s , (3.19) 
hold with u being a real number. The average ((V0)^)s is 

J2 



q^dq 1 

\Jq^ + ql 



Q = / ;L_L_ — , {aoQc « 1) , (3.21) 



where L is the system length, and the upper cutoff of the g-summation is denoted by a^^ which may be 
proportional to a~^. Luther and Peschel [20] suggested irao = a. However, care must be taken to estimate 
Q, because equation (3.20) is strongly dependent on Q. Let us consider the A = 1 (i.e., r] = 1) case. As 
discussed in §2, for A = 1, the effect of S± and Sz cancel out with each other when equation (2.6) holds. 
Here we estimate so that the equation 

{B±cose + BziS/6fcose)s = 0, (3.22) 
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reproduces equation (2.6) when A = 1, which yields 

Q = tt/o^ . (3.23) 
After some calculations, for the Haldane state, we obtain 

_..,exp(-«2s)2«||!>s. ,3.4) 

The second term of the rhs of equation (3.24) can be dropped, because ((V0)^)s is almost independent of 
B, as can be seen equations (3.20) and (3.21). Then, equation (3.15) with equation (3.24) is reduced to 

B=i«p(-<^2>i)(«. + a^.)- 
The quantity (^^)s is estimated as 

(^^)s = iE;:;^ = '^iog^' Nc«i), (3.26) 

which leads to the self-consistent gap equation 

Similar calculation can be performed for the LD state. If we replace 5±_ + {r]/2) 6z in equation (3.27) by 
\5± + {rj/2) 5z\, the gap equation becomes valid both for the Haldane state and for the LD state. Thus the 
excitation gap behaves as 

2/(4-,,) 

(3.28) 



The Haldane-LD phase boundary can be obtained from Cgap = 0, which leads to 



2(1 - I J'l) + (1 - ^ + I J'l) 1^1 + i^ Yl^jl^ j = • (3.29) 

where equations (2.3), (3.3) and (3.9) are employed. Okamoto, Nishino and Saika [10] have already obtained 
this phase boundary equation, compared it with Hida's numerical result [8] and discussed its validity. 
If the critical value of D is denoted by Dc when \J'\ is fixed, we can rewrite equation (3.28) as 



'-gap 



'\D-D,\-^, (3.30) 
2 

iy=- . 3.31 

4 - 

Thus we obtain 

m = i^B=i^ , (3.32) 

where z/H (i^ld) is the critical exponent of the excitation gap when the Haldane-LD boundary is approached 
from the Haldane (LD) phase. 

The calculation of the longitudinal spin correlation {SjSi) by the use of Hs has been already done 
by the present author [21]. Here we only write down the final expression for {SjSf), without entering into 
details; 

(5|5f ) ~ exp (^-\t^'^ , {\j - ;| oo) , (3.33) 
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where ^ is the correlation length defined by 

|£)-£>cr" ■ (3.34) 



Egap 



Thus V is also the exponent of the correlation length as expected. Comparing the result of reference [21] with 
the exact result for the XY case [22], we may rely on the exponential factor in equation (3.33), although the 
preceding power factor is not reliable. 

Hida [8] proposed the string order parameters for the Haldane phase and the LD phase. In the S 
representation, these string order parameters are defined by 

Og= hm Og(.7-0, (3.35) 

0£d= lim O£D(i-0. (3.36) 

\]-l\^oo 

with 

0^{j -l) = -4(4" eM^iS^j+i + 4" +2 + • • • + 4;-2)}4"-i) , (3.37) 

0£D(i - = -MS^j-i exp{i7r(4« + S^j^, +■■■ + , (3-38) 

a = x,y,z. (3.39) 

By the use of the identity 5| = exp(i5|)/2i valid for the spin-1/2 operators, we can rewrite equations (3.37) 
and (3.38) as 

0^(i -l) = (exp{i7r(5|,. + + • • • + S|;_i)}) , (3.40) 
0£D(i - = (exp{i7r(5|,._i + + • • • + 5|;)}) , (3.41) 

respectively. Since the slowly varying part of the z-component of the spin density is expressed as (l/27r) {d6/dx),^ 
the boson representations of 0^{j — I) and 0£j-,(j — /) are 

0^^{x-x') = Ol^ix-x') 

= {cxpmx)-9ix')]/2}) 
= (exp{i[</.(a;) - <P{x')]/2}) 

= 0^x-x'). (3.42) 

We note that the difference between the 0^{x — x') and 0£j-,(x — x') is lost due to the continuum approxi- 
mation used in the bosonization method. 

Let us calculate in the framework of the SCHA. Due to the harmonic nature of Wsi we have 

0^{x - x') = exp |-i([e(x) - 9{x')f)s^ 

= exp(^-i(<^2^s) exp|i(<^(x)<^(x'))s| . (3.43) 
The average {<j){x)(f){x'))^ is estimated as [21] 

{cj){x)cl>{x'))s = —Ko{q.[\x - x'\ + ai]) , (3.44) 

ai = , 7 = 0.5772 . . . (Euler's constant) , (3.45) 

TT 

where Kn (y) is the n-th order modified Bessel function of the second kind. We note that the parameter ai 
is determined so that equation (3.26) is reproduced when ja; — x'\ = 0. Using the asymptotic behavior of 

Ko{y) 

Ko{y) ^ J^e-« , (y 00) , (3.46) 



we obtain 

~ g?/^ ~ |£) - (3.47) 

If we denote the exponents of and 0£j-, by 2/?^ and 2/3£j-, respectively (for the factor 2, we follow Hida's 
definition [8], we see 

/3h = /3ld = = P . (3.48) 

16 — 477 

We can also calculate the asymptotic behavior of 0^(a; — x') when |a; — — > oo just on the Haldane-LD 
boundary. Using equations (3.43)-(3.45) and taking the limit — > 0, we obtain 

0'{x - x')\d=d, ~ k - a;'!""/^ , (3-49) 

where the asymptotic form 

i^o(?y)^-logy , (y^O), (3.50) 

is employed. Hida [8] also numerically calculated the system size dependence of and 0£q on the Haldane- 
LD boundary. He defined 

0^iN)\D=D^ - N-i^" , (3.51) 

Ol^mo^D^ - N-f^-^ , (3.52) 

where N is the system size. Although Hida [8] used rj^ for this exponents, we use /x to avoid confusion. From 
equation (3.49), we can expect 

Mh = t^LD = I = M • (3.53) 
This result can also obtained from the finite size scaling ansatz 

0'{\D - D^\,N) ~ N-i'fiNlD - D^n , (3.54) 

where we note that N/^ ~ N\D — -DcT- The form of equation (3.54) is chosen so that equations (3.51) and 
(3.52) are reproduced a,t D = D^. If we fix \D — Dc\ and take the limit N ^ oo, we obtain 

f{x) ~ a;'' , {x^oo) , (3.55) 
0^\D - D^\,N = oo) \D - D^l^"" , (3.56) 

from the condition that 0^{\D — Dc\,N) should be independent of N when N ^ oo. From equations (3.47) 
and (3.48) it follows that 

Hu = 2P , (3.57) 
which readily leads to equation (3.53). The relation (3.57) was also noticed by Hida [8]. 

§4. Discussion 

We have calculated several critical exponents of the Haldane-LD transition. These exponents are 
controlled by f, which varies continuously on the Haldane-LD boundary. Thus the present transition is of 
the Gaussian universality class. This fact has been already pointed out by several groups [8-10]. 

In this paper we have used the mapping of the original spin Hamiltonian onto the generalized version 
of the sine-Gordon Hamiltonian. Since several approximation have been used in the course of the mapping, 
the values of Dc and the critical exponent themselves are unreliable. We believe, however, that the present 
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theory well (lesc;ribes the qualitative feature of the Haldane-LD transition. Prom equations (3.31), (3.48) 
and (3.53), our prediction for the relations between critical exponents are 

1^ = z^D = , (4.1) 
/3h = /3ld = /? = , (4.2) 

m = mld ^ ^ ^ ^ ■ (^-^^ 

The reliability of the present theory can be examined by testing whether equations (4.2) and (4.3) hold 
or not in Hida's numerical result [8,23]. This test is summarized in Table I, where the exponents u, /3h, 
/?LD, Mh and /uld are numerical results by Hida and /?theor and /Utheor are calculated from Hida's u through 
equations (4.2) and (4.3). As can be seen, the values of /3thcor and /ithcor very well agree with /?h and /?lD: 
and /xh and /xldj respectively. Therefore we can say that our theory successfully describes the qualitative 
feature of the Haldane-LD transition. 
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Figure Captions: 

Fig.l Important parameters in Hamiltonian (2.2) on the \ J'\ — D plane. The transverse and longitu- 
dinal alternations arc mutually competing in the shadowed areas. The point S(|J'| = 1, D = 2) 
corresponds to the uniform isotropic antiferromagnetic Heisenberg model. 
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